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Motivated by experimental evidence of violations of the no-slip boundary condition for liquid flow 
in micron-scale geometries, we propose a simple, complementary experimental technique that has 
certain advantages over previous studies. Instead of relying on externally-induced flow or probe 
motion, we suggest that colloidal diffusivity near solid surfaces contains signatures of the degree of 
fluid slip exhibited on those surfaces. To investigate, we calculate the image system for point forces 
(Stokeslets) oriented perpendicular and parallel to a surface with a finite slip length, analogous 
to Blake's solution for a Stokeslet near a no-slip wall. Notably, the image system for the point 
source and perpendicular Stokeslet contain the same singularities as Blake's solution; however, each 
is distributed along a line with a magnitude that decays exponentially over the slip length. The 
image system for the parallel Stokeslet involves a larger set of fundamental singularities, whose 
magnitude does not decay exponentially from the surface. Using these image systems, we determine 
the wall-induced correction to the diffusivity of a small spherical particle located 'far' from the 
wall. We also calculate the coupled diffusivities between multiple particles near a partially-slipping 
wall. Because, in general, the diffusivity depends on 'local' wall conditions, patterned surfaces 
would allow differential measurements to be obtained within a single experimental cell, eliminating 
potential cell-to-cell variability encountered in previous experiments. In addition to motivating the 
proposed experiments, our solutions for point forces and sources near a partial-slip wall will be 
useful for boundary integral calculations in slip systems. 

I. INTRODUCTION 

Recent reports of an apparent breakdown of the no-slip boundary condition for liquid flows in small geometries 
provide an exciting and surprising opportunity to revisit one of the most fundamental questions in hydrodynamics. 
It has become textbook knowledge that, in the framework of continuum mechanics, the velocity of a viscous fluid at 
a solid boundary is equal to that of the solid. If the solid is at rest, the adjacent fluid must also be at rest. Although 
it can not be derived from first principles, decades of agreement with experiments has led to a consensus that the 
no-slip boundary condition is indeed correct for the fluid/solid boundary [1]. 

The simplest and most natural violation of the no-slip condition would involve a surface slip velocity that varies in 
proportion to the local shear rate, written in the case of a flat surface as 

u,|=A-^i u _l = 0. (1) 

This condition naturally introduces a new length scale, A, called the slip length. In gases, the slip length is related to 
the mean free path, A / , where non-continuum effects become important [2] . The analogous picture does not appear 
to hold for liquids, however. Liquid molecules are in constant collision, and any analogous mean free path would 
be of molecular order, significantly smaller than recent experiments suggest (discussed below). Regardless of the 
discrepancy between physical origins of apparent slip in liquids and gases, the effects of slip are expected to become 
important when the experimental length scale h is of the same order as the slip length A. Therefore, by analogy with 
gaseous slip flows, we will use an effective Knudsen number, 

Kn=£ (2) 

to describe flows near a partial-slip surface. Obviously, we expect slip effects to play a significant role when Kn > O (1). 

In the past half-century, the no-slip condition saw only occasional challenges [3, 4]. More recently, however, various 
experimental systems have probed liquid flows on small enough length scales h that Kn may no longer be small, 
allowing a more thorough and sustained re-investigation of the no-slip boundary condition. These experimental 
techniques differ in the way that flow is created and slip is measured, and fall into five primary categories, reviewed 
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in [5]. (i) One can measure the relationship between flow rate and pressure drop in capillaries or microchanncls 
[6, 7], which depends on Kn. (ii) One can measure the force required for squeeze flows in long and narrow geometries 
such as are found in the surface force apparatus (SFA) or atomic force microscope (AFM) [8-19]. Slip reduces the 
viscous resistance, and quasi-steady probe motion is assumed, (iii) One can measure pressure-driven velocity profiles 
in a capillary or microchannel using small particles as passive tracers [20, 21]. Tracers that are sufficiently small, 
uncharged, and far from the wall should faithfully reproduce the fluid velocity, although the high diffusivity of small 
probes requires averaging techniques, (iv) One can measure an externally-driven flow near a wall using fluorescence 
correlation spectroscopy with labelled molecular tracers [22]. (v) One can use near field laser velocimetry, wherein 
evanescent optical waves (exponentially localized to a small region near a wall) are used to measure the velocity of 
photobleached molecular probes in an externally-driven flow [23, 24]. 

There are now many published reports of apparent violations of the no-slip condition, both experimental [6-27] and 
theoretical [28-30]. Apparent slip has been measured over surfaces that are completely wetting [9, 10, 23], partially 
wetting [12, 13], and non-wetting [7, 8, 11, 13, 20, 30]. Roughness has been predicted and measured to decrease 
slip [14, 23, 31, 32], although in some cases roughness appears to increase slip [10, 33]. In some measurements and 
simulations, the slip length appears to be independent of shear rate [8, 11, 23, 30], whereas in others it depends upon 
shear rate [7, 9, 10, 12, 13, 29]. Moreover, apparent slip lengths ranging from nanometers [9] to microns [13] have 
been been reported. It is thus reasonable to conclude that no consensus has been reached concerning the existence 
and physical origin of fluid/solid slip, and the physical factors that influence it. 

The large variability in the published results could be due in part to the variety of experimental techniques employed. 
After all, physical mechanisms other than liquid/solid slip can resemble apparent slip in experiments [34-38], and could 
lead to incorrect conclusions as to the nature of the actual solid/liquid interface. Different experimental techniques 
are susceptible to these effects to different degrees. Additionally, all require an externally-forced flow or motion, which 
introduces an additional source of experimental uncertainty. Furthermore, multiple experimental cells are typically 
required to probe different solid/liquid surfaces. Finally, many experiments involve averaging - over the length of a 
capillary, the area of an SFA, or the diffusive motion of tracers. 

In this paper, we propose a complementary technique to probe the nature of the liquid/solid boundary that is 
largely immune to the issues raised above. The idea is to measure the influence of the wall on the Brownian motion of 
suspended tracers. A spherical particle of radius a, far from the wall, diffuses with a bulk diffusivity Dq = ksT/GiTfia. 
When a wall is located a distance h from the particle, particle diffusivity is affected in a manner that depends on the 
nature of the surface. A no-slip wall (Kn = 0) gives corrections to the perpendicular (D ± ) and parallel (£)") diffusion 
coefficients, 



with errors of order 0(a 3 /h 3 ) [39]. If, however, the surface is perfectly slipping (i.e. sustains no shear stress, or 
Kn = oo), the particle diffusivities are given by 



Finite values of the slip length (or Kn) should interpolate between these two limits. We note in particular that the 
parallel diffusivity goes from being wall hindered for Kn C 1 to enhanced for Kn 3> 1. Naturally, a knowledge of 
the relation between slip length and diffusivity would allow the slip length of a solid/fluid interface to be inferred 
from measured diffusivity of nearby particles. No external flow is required, and walls with patterned wettability allow 
various surfaces to be probed within a single experimental cell, which would allow differential measurements that are 
free of the uncertainties due to cell-to-cell variability. This builds upon an idea that was first pursued by Almeras et 
aJ. [40] , who characterized the influence of wettability and slip on the parallel diffusion coefficient of a small particle 
between two walls, as involved in molecular diffusion under confinement. 

In this work, we calculate fundamental solutions for Stokes flows near a single partial-slip wall. Our results give 
an explicit relationship between solid/liquid slip and colloidal diffusivity, as well as expressions for the flow fields 
themselves. In addition to aiding in intuition for partial-slip systems, the flow fields we calculate will be useful for 
boundary integral calculations in partial-slip systems. Additionally we explore the feasibility of measuring the effect 
of a partial-slip wall upon colloidal diffusivity as a means of measuring the wall slip itself. Recent years have seen 
precise experimental measurements of colloidal diffusivity near walls and/or other colloids. Corrections of order a/h 
can be accurately measured, and excellent agreement has been found with theory [41-43]. 

This technique has various advantages. First, it does not require an external flow and therefore alleviates the ex- 
perimental difficulties associated with precise flow manipulation. As a consequence, the experiment can be performed 
in a closed cell, and thus avoid contamination by impurities. Consequently, the liquid can be degassed or put under 
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FIG. 1: First image of the Stokeslet; (a): Stokeslet (F, h) perpendicular to the surface; the first image is the Stokeslet (-F, —h); 
(b): Stokeslet (F,h) parallel to the surface; the first image is the Stokeslet (F, —h). Note that these would be the complete 
image systems if the surface was perfectly slipping (Kn = oo). 

variable pressure to probe the influence of adsorbed nanobubbles, as discussed below. Second, our method does not 
average over different experiments, sample volume or apparatus size but instead makes use of a single colloidal probe. 
Third, multiple solid/liquid interfaces can be probed within a single experimental cell by using deliberately pat- 
terned surfaces. This would allow differential measurements to be performed, and possibly to track surface-attached 
nanobubbles [44-46]. 

The paper is organized as follows. In §11, we consider the effects of a partial-slip wall on the two main fundamental 
singularities of Stokes flow - the point force (Stokeslet) and the point source. In a manner analogous to Blake's image 
system for a no-slip wall [47], we interpret the wall's contribution in terms of a series of image singularities. This is the 
central result of our work. In §111, we use this solution to provide an analytical formula for the influence of a partial- 
slip wall on the diffusivity of a small spherical particle. In §IV, we consider the coupled mobilities/diffusivities of two 
small particles, and propose alternate experimental tests for slip to complement those in §111. As the calculations 
themselves are somewhat laborious, we relegate the details to appendices, and save the main body of the text for key 
results and discussion. 



II. IMAGE SYSTEMS NEAR A PARTIAL SLIP SURFACE 

Analogous to point charges and point masses in electrostatics and gravitation, flow fields associated with fundamen- 
tal singularities are useful in treating Stokes flows [48] . There are two families of fundamental singularities in Stokes 
flows - (i) the point source/sink and their derivatives, which correspond to irrotational potential flow and are entirely 
analogous to electrostatic fields, and (ii) the point force (Stokeslet) and its derivatives, whose flow fields are viscous 
and rotational. These fundamental singularities aid in intuition for viscous flows, and in providing approximate and 
asymptotic solutions. Furthermore, they form the basis for boundary integral techniques in Stokes flow calculations 
in more complicated geometries, where flow and pressure fields are computed by solving for surface distributions of 
fundamental singularities [48]. 

In a classic paper, Blake [47] interpreted the flow field due to a Stokeslet near a no-slip surface in terms of a 
system of image singularities, located on the opposite side of the wall. Blake's image system consists of an equal but 
opposite Stokeslet, a Stokeslet dipole (i.e., force-dipole) and a source dipole (potential dipole). A perfectly slipping 
surface (Kn = oo), has a simpler image system: A single Stokeslet of equal magnitude and symmetric direction, as 
in Fig. 1. The image system for a partial-slip surface (0 < Kn < oo) is more complicated, and is the subject of the 
following analysis. In what follows, we present the complete image system for a Stokeslet and a point source near a 
planar partial-slip boundary. Higher order (multipolar) singularities can be derived from these two by differentiation, 
although subtleties exist (discussed below). Because our calculation is analogous to Blake's, but algebraically more 
involved, we save the details for Appendices A, B, C and D. 

A. Set up and boundary conditions 

We choose the x — y plane to lie along the solid wall, with the z coordinate directed perpendicular to the surface, 
and consider a Stokeslet of strength F located at (x, y, z) = (0, 0, h). The velocity field, u, satisfies the incompressible 
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Stokes equations 

fiV 2 u = Vp, V • u = 0, (5) 

subject to partial slip boundary conditions (Eq. 1). We decompose u into three components, 

u = U + V + w, (6) 

where U is the flow field due to the Stokeslet itself, V is the flow field due to the primary image Stokeslet of strength 
F located at (x, y, z) — (0, 0, —h) (Fig. 1), and w is an as yet unknown velocity field that solves Eq. (5). The Green's 
function for the Stokeslet is given by [47, 48] 



8ir[i \r 

Since the wall breaks the isotropy of the particle mobility, we consider the perpendicular ( ± ) and parallel (H) 
Stokeslets separately, giving velocity fields 

= FGf(r) (8a) 
U" = FGf(r), (8b) 

with r = (x,y,z — h). Here we have introduced the notation 

Gf(r) = G s (r-r )-e z , (9) 

where ro is the location of the Stokeslet, r is the observation point, and the subscript indicates the direction of the 
force. In what follows, we will also use the following notation for higher-order singularities 



d 



dx r =0 



d 



G °» = — { GS ( r - r o) • e z } ee -_ |G s (r) • e,} (10a) 



Gf - (r) = dSy- - ro) ' e '> S ^ ^ S(r) ' e '> • (1 ° b) 



r =0 



dxdy 



Here G SD represents a Stokeslet doublet, G s< ^ a Stokeslet quadrupole, and so on. Note that derivatives are taken with 
respect to the singularity location (rather than the observation point), one derivative for each coordinate following 
the semicolon. 

For the primary image Stokeslets located at z = —h, we pick F as in Fig. 1 to enforce the no-flux condition at the 
wall, giving velocity fields 

V 1 - = -FGf(f) (11a) 
V" = FGf(r) (lib) 

where r = (x, y, z + h). 

Enforcing the partial-slip boundary condition (Eq. 1) at the surface z = imposes boundary conditions on w, 

l-^)w,| =-2U,|, w 3 =0. (12) 

We use Fourier transforms in the x- and y-directions and the general solution to Stokes equations, as given in Appendix 
A. The amplitudes of the Fourier components are determined by enforcing the slip boundary condition (12); for details 
see Appendix B for the perpendicular case, and Appendix C for the parallel case. These solutions can be used to 
calculate the change in the particle mobility, and therefore diffusivity. 

An interesting result is that the Fourier coefficients for both cases can be related directly to the coefficients in the 
no-slip series, which allows the partial-slip solution to be expressed in terms of weighted integrals of the no-slip image 
systems. This results in a clear physical interpretation of the partial-slip image system in terms of weighted integrals 
of fundamental singularities. 



FIG. 2: Streamlines for a Stokeslet oriented perpendicular to a partial-slip wall, with (a) Kn = (Blake's solution, no-slip), (b) 
Kn = 1, (c) Kn = oo (perfect slip). The streamlines are displayed in the plane which includes the Stokeslet and is perpendicular 
to the nearby surface. 



B. Stokeslet perpendicular to slip surface 



As shown in Appendix B, the total velocity field for a Stokeslet oriented perpendicular to a partial-slip wall can be 
expressed as = U 1 - + V 1 - + w^, where XJ 1 - and V 1 - are given by Eqs. (8a) and (11a), and w 1 - is given by 



Fh I' 00 

w>, A) = — e- s / 2X [(h + s) Gf - Gf£ ](r + (h + s)e z )ds, 



(13) 



where G D is a potential (source) dipole, defined by 

1 d 



g d = 1 d ( r - r 



8ir[i dzo \ |r — r | 3 



r =0 



87r/x dz 
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Analogous formulae and notation for source quadrupoles G® and so on follow in a straightforward fashion. This 
image system therefore represents a weighted line integral of source dipoles and Stokeslet dipoles, whose magnitude 
decays exponentially with distance (scaled by the slip length). In fact, Eq. (13) represents a line integral of Blake's 
image system for no-slip walls [49]. The streamlines for the complete image system are displayed in Fig. 2. 
Scaling the integration variable s with Am gives a form for the image system, 



/>oo 

w ± (r,X) = Fh / e- u ' 2 [{h + Au)Gf - Gf.f](r+ (h + Xu)e z )du, 
Jo 



(15) 



which is amenable to asymptotic analysis. The basis for the analysis that follows is that the integrand is only 
appreciable when u < O (1), and is exponentially small otherwise. 

Since z and h play the same role in the argument of the two singularities in Eq. (15), we can consider the limit 
where A <C (z + h) in the singularities in Eq. (15). A Taylor expansion of the term in brackets gives 

w x (r) « 2Fh[(h + 2A)Gf - Gf ; f ](r) + 4Fh\[-(h + 4A)G« + G s z % ](f). (16) 

In the limit when the particle is farther from the wall than the slip length (A -C h, or Kn <C 1), Eq. (16) becomes 

w^r, A « h) = 2Fh[hGf - Gf.f ](f) + 4Fh\[G? - hGg + Gf.« ](f). (17) 

The first term is Blake's solution for a no-slip wall, and the second term represents an O (Kn) correction to the image 
system due to slip. Furthermore, in the limit where Eq. (16) results in 

w-^r) 4F/iAGf (f), (18) 



which differs significantly from Blake's solution. 



FIG. 3: Streamlines for a Stokeslet oriented parallel to a partial-slip wall, with (a) Kn = (Blake's solution, no-slip), (b) Kn = 1, 
(c) Kn = oo (perfect slip) . The streamlines are displayed in the plane which includes the Stokeslet and is perpendicular to the 
nearby surface. 



C. Stokeslet parallel to slip surface 



As shown in Appendix C, the total velocity field for a Stokeslet oriented parallel to a partial-slip wall can be 
expressed as = U" + V" + w", where U" and are given by Eqs. (8b) and (lib), and w" is given by 



f r°° 

w«(r,A) = J J o l-G s x +hG s z °-h 2 G°](r+(h + S )e z )e-°/ 2X ds 



- AFX / G™{T+{h + s)e z ) 
Jo 



e s/2X _ 1 



ds 



(19) 



+ 4FA j [Gf - hGg}(r + (h + s)e z ) [l - (l + ±) e 
where G^.® is a rotlet dipole, given by (see details in Appendix C) 

x;yy>- 



s/2\ 



ds, 



pi?D _ t(c SD r<SD 

z;y 2 y> x v 



(20) 



The solution, Eq. (19), can again be interpreted as a line integral of the fundamental singularities above, but their 
weight does not systematically decay exponentially away from the image location. These are shown in Fig. 3. Note 
that, as A — > oo, each term in Eq. (19) goes to zero due to the rapid spatial decay of the singularities (G D ~ 1/r 3 , 
G RD ~ 1/r 3 , G Q <~ 1/r 4 ) and due to the vanishing value of the weights in Eq. (19) as s — > 0. 

Re-scaling s by A in Eq. gives 



/•OO 

w" (r, A) = F / [-Gf + hG s z D x - h 2 G° ](r + (h + Xu)e z ) e^ 2 du 
Jo 



- 4FA / G™(r+(h + \u)e z ) 



-u/2 



du 



(21) 



+ 4FA 2 ^°°[Gf - hGg](r +(h + \u)e z ) [l - (l + |) e""/ 2 ' 



du, 



which leads to asymptotic formulae for the image system using an expansion of the terms in Eq. (21). In the limit 
Kn <C 1, the image system is found to be given by 



w« (r, A) « 2F[-G S X + hG s z D x - h 2 G°](f) + AXF[G S X D Z - hG s z % + h 2 Gg](r), 



(22) 



which is Blake's solution for a no-slip surface plus an O (Kn) correction. Note that the rotlet dipole in Eq. (21), G z .y, 
will only appear in asymptotic formulae for the image system at order O (Kn 2 ). 
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D. Image system for a point source 

We finally consider in this section the image system for a point source, which we will denote by In this case, 
the velocity field and its first image are given by 

u« = j£-4> ( 23a ) 

V» - (23b) 

where 47r<5 can be interpreted as the source flow rate. Here again, we decompose the velocity field u^^ = U^' + V^' + 
w( ) and need to solve for w( ). The boundary conditions for w^^ are exactly proportional to those for w 1 - (Eq. Bl); 
correspondingly the two solutions are proportional as well, giving 



c "i 2Q i 2Q 

w w = -w = — 

Fh A 



/'OC 

/ e-°/™[(h + s)G? -G s z °}(r+(h + s)e z )ds. (24) 
Jo 



E. Higher-order singularities 

Finally, we note that higher-order singularities may be derived from the Stokeslet and point source image systems 
presented above, by taking derivatives with respect to the singularity location. Subtleties do exist, however, so care 
should be taken [49]. Derivatives along the plane of the wall can be computed in a straightforward fashion, but 
derivatives perpendicular to the wall are more subtle, because the amplitude of the image singularities depends on 
the distance from the wall. In general, the correct image system will always be obtained if the derivative is calculated 
by taking the limit of two such opposing singularities (since each solution obeys the correct boundary condition on 
the wall in the first place). That is, if u sm (h) is the fundamental singularity located a distance h from the wall, and 
u l (h) is the proper image, then the image system for a perpendicular dipole of u sm (/i) can be found by taking the 
limit 

U D = Um / u si » + u'(fe) - u^(h - e) - u'(fe - e) \ _ (25) 



III. INFLUENCE OF SLIP ON BROWNIAN MOTION 
A. Diffusion of a spherical particle 

We now turn to a specific application of the above calculation: the diffusivity D of a solid spherical particle of 
radius a near a partial-slip wall. Using the Stokes-Einstein relation, the diffusivity is directly proportional to the 
particle mobility b via D = fcgTb; thus the (deterministic) calculation of mobility yields the diffusivity. 

Particle mobilities are defined as the velocity response to a force F acting on the particle. In the absence of solid 
boundaries, and if the no-slip boundary condition is satisfied on the particle surface, the velocity field established 
around the particle is given by 



/ 1 rr \ „ a 2 ( 1 3rr > 

u = o - + -5 ' F + — r ' F ' 26 

8irn \r r 3 ) 24tt^ V ^ 3 r ' 



leading to the (isotropic) Stokes mobility 



' =>0 Gir^a' 



The first term in Eq. (26), which decays like 1/r, corresponds to a Stokeslet, whereas the second term (source dipole) 
is necessary to satisfy the no-slip boundary condition on the surface of the particle, and decays like 1/r 3 . 

The presence of a nearby surface (a distance h from the particle) modifies the flow field around the particle, and 
hence its mobility, which is now an anisotropic tensor b. To account for such effects, one can employ an approximate 
and iterative technique known as the method of reflections [39]. When the particle is 'far' from the wall (a -C h), it 
sets up a flow field that is locally appears like that around a particle in an infinite fluid. This flow field, however, 
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violates the boundary conditions at the wall, and so an additional field ('reflection') is introduced to correct the 
boundary conditions at the wall. This first reflection, however, violates the boundary conditions at the particle 
surface, necessitating a second reflection, and so on. 

Thus in our approximation, the particle travels through its local fluid environment with Stokes mobility (27), and 
is advected by the image flow u w via Faxen's law 

n 2 

u adv = u w + =-V 2 u w . (28) 
6 

Thus the leading order correction (in a) to the mobility is given by the wall-induced flow u w evaluated at the particle 
location. Furthermore, to obtain the O (a) component of u w , only the image system for the Stokeslet in (26) is 
required. Errors to this approach are of order a 3 /h 3 , since the source dipole flow in (26) and the Laplacian in (28) 
are smaller by a 2 . Note also that, since only the image system for the Stokeslet is required for the leading-order wall 
correction to the mobility, the result is insensitive to the boundary conditions (i.e. size, shape, or slip) on the particle 
itself. 



B. The effect of slip upon single-particle diffusivity 

The components of the velocities at the position of the Stokeslets and in the same direction as the applied force, 
denoted generically (ui, 112,113), are given by 

uj = ^- + v 3 ± (0,0,h) + wi(0,0,h), (29) 
+ V} l {0 7 0,h) + w l l{0 7 0,h). (30) 



6717x0 

The first terms (U~ L ,U") represent the bulk (Stokes) mobility; the second terms (V^, V") represent advection with 
the flow field established by the primary image Stokeslets, and the third terms (w 1 ^') reflect advection with the 
higher-order image field, as detailed in Appendices B and C. The flow velocities from the primary image Stokeslet, 
evaluated at the particle position, are 

^<°'°'*> = -8^' *°'*> = I^' (31) 

and inverting the Fourier transforms from Appendices B and C reveals the contribution from higher-order singularities 
to be 

wi(0,0,h) = -^J(Kn), (32a) 
wl(0,0,h) = -J^J(Kn) (32b) 



where the functions X and J arc defined by 



x 2 _ 2 _ Kn(Kn - 1) + e^ Kll T(0, 1/Kn) 



I<K °» " L Y+2xKi ' " — 8Ktf ' ■ < 33i " 

V ' J (l + Knr)(l + 2Knr) V ; 

-Kn(3Kn + 1) + e 1 / Kn (l + 2Kn) 2 r(0, 1/Kn) e 2 / Kl T(0, 2/Kn) 



8Kn 3 Kn 

where T is the incomplete Gamma function [50] 



/>oo 
J X 



T(a,x) = / t a -'e- l At. (34) 
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FIG. 4: Variations of X and J with the Knudsen number Kn. 



The results of Eqs. (29) and (30) yield the desired relation between wall slip and nearby colloidal diffusivities, given 
to leading order in a/h by 

£-l-£(l + *r(K=)) + o(£). (35.) 

£- 1 + £ (1 - 17 (Kn)> + £>(£)■ (35b) 

The functions I and J" are illustrated in Fig. 4. It is significant to note that appreciable variation in parallel diffusivity 
(J) occurs over about three decades in Kn, giving a fairly wide range of experimental conditions under which one 
might hope to measure A. 
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C. Asymptotic limits 

When the slip length is small (A <C h), we can use a Taylor expansion of Eqs. (35a) and (35b) and obtain the 
asymptotic results 

£i = 1 _g (l _ Kn + 0( Kn>)) +0 (£), (36a) 



£--&(i-k. + ^-) )+0 (£). 



(36b) 



The results for a no-slip surface, Eq. (3), are therefore recovered as the slip length vanishes. The 'slightly slipping' 
result can also be obtained simply from the no-slip solution, using a reciprocal theorem for Stokes flow, as shown in 
Appendix E. 

When the slip length is large A » h, on the other hand, Eqs. (35a) and (35b) can be approximated by 

The results for a no-shear surface, Eq. (4), are recovered when the slip length diverges, albeit slowly (note the 
logarithmic dependence in the parallel case). 
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IV. COUPLED DIFFUSION OF TWO PARTICLES NEAR A PARTIAL SLIP SURFACE 



We extend in this section the idea proposed in §111 to the coupled diffusivity of two colloids. We consider two 
spherical particles, radius a, located at distance d from each other along the z-axis and at a distance h above the 
slipping surface. 

The diffusivity tensor, D, for a general N-particle system is given by [51] 



D = k B Tb 



(38) 



where b is the N-particle mobility tensor. In the case of two particles, these mobilities are the tensors relating the 
forces, Fi and F 2 , acting on each particle and their velocities, Ui and u 2 , as 



Ui — bn • Fi + bi2 • F2, 
u 2 = b 2 i • Fi + b 22 • F 2 . 



(39a) 
(39b) 



The values of the tensors bn = b 22 reflect the self-diffusion of each particle, and the influence of slip on their values 
was presented in §111. We calculate below the influence of slip on the coupling tensors, b^ 2 = b 2 i. Note that unlike 
bn and b 22 , the coupling tensors bi 2 and b 2 i do not depend on the particle size a, but rather on the value of the 
new length scale, d. 

Since there is symmetry between particles 1 and 2, we have hf 2 = b 2 i, and it is therefore sufficient to calculate 
only b 2 i. Furthermore, since the particles are aligned along the x-direction, the tensor has only five non-zero entries, 



h 



b 

1/22/1 





X2Z1 




(40) 



By symmetry, b X2Zl = —b Z2Xl , leaving just four independent components in the coupling mobility/diffusivity. 

We adopt the same method as in §111 and consider the limit a <C h, so that we replace the full velocity field around 
the particle by a Stokeslet. In the case of no-slip, the coupled mobilities (influence of the point force and its image 
system) are given by 



b 21 (Kn = 0) = 



2 1- 



V 



(i + B/2 


3^3/2 



(i + £) 5/2 

whereas when the slip length is infinite, they are given by 



b 2i (Kn = 00) = 



(l + £)3/2 



£1/2 



V (i + 3/2 



with 







1 



(l + £)3/2 




_3£3/2 
(l + 5/2 





\ 



1 - 







1 



(1+0 1/2 





(1+^)3/2 





1 - 



1 + f j + 3£ 2 

(i + 5/2 ) 



£1/2 \ 



1 + 2g 

(i + 6 3/2 / 



(41) 



(42) 



4^2 
?_ d? ' 



(43) 



These results are correct to order O (a 3 //i 3 ) and O (a 3 /d 3 ). The qualitative difference between the no-slip and perfect 
slip formulae concern their spatial decay in the limit £ — ► (h -C d). In that case, all components of the mobility 
tensor decay faster in the case of no-slip than in the case of perfect slip. 

We present below the calculation for b Z2Zl , as it is the configuration where calculations are the easiest. Calculations 
for other components of the mobility tensor are derived in Appendix F, with results summarized below as well. 

In order to evaluate b Z2Zl , we consider a unit vertical force applied to particle 1, and determine the vertical velocity 
at the position of particle 2. Two factors contribute to the value of b Z2Zl : the direct influence of the Stokeslet flow 
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field from particle 1 and the influence of the image system for this Stokeslct below the slipping surface. Given the 
decomposition assumed in Eq. (6), the mobility is given by 

b Z2Zl {Kn) = b Z2Zl (Kn = oo) + iof(de x + he z , A). (44) 

Using Eq. (42) and and evaluating the integral in Eq. (13) leads to 



b z 2Zl (Kn) 



1 - , 1 + 2 L + U J- f e -VKn 1-2^(1 + ! ; du 



,2 



(l+0 3/2 2 \KnJ (l + e(l + «) 2 ) 5/2 



(45) 



For a fixed value of Kn, the asymptotic behavior of the integral in parenthesis in Eq. (45) in the limit £ — > is given 

by 



-f 

Kni 



u/Kn JTT§T^W du = 1 " ^ (1 + 2Kn + 2Kn2) * + ^ (46) 



Consequently, for any value of the slip length, there exists a range of £, 

£ « min (l, -^j , (47) 

which corresponds to the far-field limit 

d 2 > max(/i 2 ,A 2 ), (48) 

for which the integral in Eq. (46) goes to one asymptotically. It follows that, evaluating Eq.(45), the mobility always 
decays asymptotically as O (£ 2 //z(i) ~ 1/d 5 ; this is the same power law as the no-slip case. The case of perfect slip 
is therefore a singular limit: as Kn — > oo, the range of £ for which this asymptotic behavior is valid, £ <C 1/Kn 2 , 
shrinks to zero, resulting in an asymptotic behavior of the perfect-slip diffusivities qualitatively different from that of 
any other partially slipping surface. We can then use these results to obtain, for a given slip length, the asymptotic 
behavior of the mobility as £ — > 0. Substituting the result of Eq. (46) in Eq. (45), we obtain 

b Z2Zl (Kn) = - g^(l + 4Kn + 4Kn 2 ) [1 + 0(0] . (49) 

As is obvious in Eq. (49), the numerator is a function of Kn, and as consequence, a measure of the behavior of the 
spatial decay of the coupled diffusivities of the two particles allows, in principle, to infer the value of the slip length. 
The other components of the mobility tensor b 2 i are calculated in Appendix F for small values of £. They are given 

by 

WKn) = _^L(l + 2Kn + Kn 2 ) [1 + 0(0], (50) 
Kyi ( Kn ) = (1 + 4Kn + 12Kn 2 + 16Kn 3 + 8Kn 4 ) [1 + O (£)] , (51) 

b X2Zl {Kn) = -^^(l + 3Kn + 2Kn 2 ) [1 + 0(0]- (52) 

and similarly, their measurement would allow an estimation of the slip length on the surface; Eqs. (50) and (51) are 
valid when d ^> max(/i, A) and Eq. (52) when d 2 ^> max(/i 2 , A 2 ). Note that the largest leading-order influence of a 
non-zero slip length are obtained for the components b Z2Zl and b V2Vl of the coupled matrix (behavior w 1 + 4Kn for 
small Kn), although these components decay most quickly with £. The most slowly-decaying coupling mobility is 

bx'yXi • 



V. DISCUSSION AND CONCLUSIONS 



Having presented the fundamental singularities for Stokes flow near a partial-slip planar wall and explored the 
consequences for colloidal diffusion, we now turn to examine issues relevant to experimental studies of slip. 

Before we begin, it is significant to note that there are several length scales inherent in the systems we have been 
considering: the colloidal radius a, distance from the wall h, slip length A, and (for multi-particle systems) the 
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distance d between the particles. The theory presented above concerns the fluidic response to a point force, and 
thus the results we presented are valid for particles that are 'far' from the wall, so that a <C h. (Note, however, 
that excellent agreement between theory and experiment was obtained for the coupled colloidal diffusion even for 
systems with a/h ~ 1/3 [41].) This is significant because the transition between 'no-slip' and 'perfect slip' occurs, 
not surprisingly, around Kn ~ O (1), or h ~ A. Thus to experimentally observe the transition between the two slip 
regimes (and thus measure A convincingly), a tracer is required that is of the same order as the slip length itself. 

Different slip lengths could be probed with different experimental techniques. Several experiments [4, 13, 20, 22, 25] 
have reported slip lengths of order microns, which would allow the use of micron-sized colloids. Optical tweezers 
can trap colloids of this size, and thus allow their repeatable and precise three-dimensional placement. Repeatedly 
trapping and releasing single or multiple colloids has proven an excellent method for measuring spatially-varying 
single- or multi-particle diffusivities [41-43], and would thus be naturally adaptable to probe the slip properties of 
walls as described above. Such studies have been developed using video microscopy, which most typically would 
measure motion parallel to the wall, and project out perpendicular motion. Again, this requires a probe that is large 
enough to be trapped by optical tweezers, which places a lower limit on the slip length that could practically be 
measured, although a null measurement using tweezers would be useful in putting an upper bound on the slip length 
of a surface. 

Other slip experiments [3, 6-8, 10-12, 14-19, 21-23, 26, 27] report shorter (10 — 100 nm) slip lengths, which would 
require smaller (10 — 100 nm) tracers. For visualization, such tracers should presumably be fluorescent, such as quan- 
tum dots (see, e.g., [22]) These small colloids are more difficult, if not impossible, to hold with optical tweezers, which 
rules out their precise manipulation and placement. Instead, techniques such as total internal reflection microscopy 
allow accurate three-dimensional measurements of their positions, and could thus be used to probe surfaces with 
10 — 100 nm slip lengths. Such techniques have the additional advantage of precise, three-dimensional measurements, 
and could thus be used to probe both the perpendicular and parallel diffusivities. Additional issues also arise with 
such small tracers: obviously, the diffusive motion is significantly higher, so correspondingly particle motion must be 
resolved on faster time scales. Furthermore, colloid/wall interactions (electrostatic, van der Waals, and so on [52]) 
can become significant at these shorter length scales, and must be treated properly in data analysis. 

We now discuss issues specific to single-particle diffusivities. From Eqs. (3) and (4), it is evident that variations in 
Kn give rise to a change in single-particle diffusivities of at most 

AD 1 _ 3a AD" _ 15a 

to leading order in a/h. The effect on the diffusivity parallel to the wall is larger than for the perpendicular diffusivity, 
and furthermore changes sign: parallel diffusivity is hampered by a no-slip wall, but enhanced by a significantly slipping 
wall, with a crossover occuring at Kn « 5.45, where j7(Kn) = 1/2. Furthermore, from Fig. 4, the parallel diffusivity 
changes over about three decades in Kn, providing a further advantage to the parallel mode of measurement. In 
principle, significantly smaller particles (a <C h ~ A) could be used in such measurements; however, the correction to 
the self-diffusivity is of order a/h, and thus the smaller the probe, the smaller the effect to be measured. 

Multi-particle diffusion, on the other hand, affords a greater variety of measurements. In [43] and [41], pair diffusiv- 
ities were measured in terms of centcr-of-mass and relative motion variables, wherein corrections due to hydrodynamic 
interactions are smaller than the bulk 'Stokes' diffusivity by factors of order a. Rather than measuring these particular 
modes, we suggest simply measuring the cross-correlation between two probes: 

~{A ai &/h) = k B Tb aifh , (54) 

where a and (3 can take the values {ir, y, z} and 6 represents the coupling mobility. Notably the coupling mobilities 
depend on neither the shape nor the size of the colloids themselves. (This feature plays a significant role in so-called 
two-point microrheology [53].) From Eqs. (49)-(52), one can see that the largest effect of slip upon cross-correlated 
diffusive motion occurs for the zz and yy modes, although the xx has the slowest spatial decay and may be easiest to 
measure (see Fig. 5). 

Another attractive feature of the proposed experimental system is the ease of performing multiple experiments 
within the same experimental cell. Surfaces with patterned properties could be used to probe surfaces with different 
putative slip lengths, and differential measurements could be used to remove the uncertainties associated with cell- 
to-cell variability. 

Although the calculations presented here are valid only for particles 'far' from the wall, the experiments proposed 
here need not be performed in this limit. In fact, the results we present here could assist in boundary- integral studies 
to obtain the various mobilities/diffusivities for systems with colloidal spheres 'near' partially-slipping walls. 

Finally, although the idea proposed here is concerned primarily with passive microrheology, the change in the 
particles mobilities could in theory also be measured using active microrheology, i.e., measuring the direct relationship 
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d/h 

FIG. 5: Variation of the coupled mobility b X2X1 (non-dimensionalized by 1/4-K^ih) with d/h, for four values of the Knudsen 
number: Kn = 0.1 (circles and solid line), Kn = 1 (squares and solid line), Kn = 5 (triangles and solid line) and Kn = oo 
(dashed-dotted line). Insert: same data but in log-log scale, with included the asymptotic behaviors for each value of Kn < oo 
as given in Eq. (50) (dashed lines). 

between particle motion and a known applied force, whether that force were applied to the particle itself (in which 
case the self-mobility would be measured), or to an adjacent particle (in which case the coupling mobility would be 
measured) . 
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APPENDIX A: GENERAL SOLUTION OF THE STOKES EQUATIONS NEAR A PLANAR BOUNDARY 

We present in this section the general solution of Stokes equation near a solid boundary, using the notations 
introduced in §11 A, for a general velocity field u. We introduce the Fourier transform in the directions parallel to the 
plane, for each of the velocity components 

u 3 {k 1 ,k 2 ,z)=F[u 3 ] = ±- f fu j (x ) y,z)e ikix+ik ^dxdy, (j = 1,2,3) (Al) 



and its inverse 



2tt 

The general solution to Eq. (5) is given by 



u t (x,y,z)=?- i lu j ] = -l- I fn t (k 1 ,k2,z)e-' k "-' k "dkidk 2 . (A2) 



(— 

(— ) 



{A a +ik a Bz}e~ kz , a = {1,2} (A3a) 



u 3 = ^—]{A s + kBz}e-* a (A3b) 
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where k = {k\ + k 2 ) 1 / 2 and a can take the values {1,2}, and where we have taken the magnitude of the force to 
be F = 1 to simplify the notations (as the equations are linear with respect to F, this can be done without loss of 
generality). The four dimensionless constants (Ai, A 2 , A 3 ,B) are linked through the continuity equation, giving 

+ k 2 A 2 ) = k(B - A 3 ). (A4) 

The remaining three constants are found by applying the slip boundary conditions, Eq. (12), on the surface, which 
determines the velocity field uniquely. 

We note for future use that the Fourier transform of the velocity field (A3) can be divided into two components: 

u = Ua + zii hl (A5) 

giving corresponding real-space velocity fields 

u = u a + zu b , (A6) 

where u a = ,F _1 [u a ], and u b = JF _1 [u h ]. 



APPENDIX B: SOLUTION FOR PERPENDICULAR STOKESLET 



The boundary conditions, Eq. (12), for the flow w 1 - due to a Stokeslct perpendicular to the surface are thus given 
by the no-flux condition, = 0, and the slip condition 



;f - A 



~dz~ 



(x,y,0) = -2U£ = ±(^,]L), 
Anfi \r 3 h r\) 



where r\ = x 2 + y 2 + h 2 and a = {1, 2}. In Fourier space, these conditions become = and 



- A 



dw 



7,-1-1 



8z 



(k x ,k y ,0) - i^-jr? kh - 



Using the general formalism of Eq. (A3) together with the boundary conditions Eq. (B2), we obtain 

2ik-i 



At 
B 1 - 



-kh 



k(l + 2\k) 

2ih 



k(l + 2\k) 
0, 

-2 



2 e~ k \ 



kh 



l + 2\k 



Writing these constants as 



Ai(k,z,X) = 
B^(k,z,X) = 



^(k,z,0) 
l + 2Afc ' 

l + 2Xk ' 



(Bl) 

(B2) 

(B3a) 

(B3b) 

(B3c) 
(B3d) 

(B4a) 
(B4b) 



where Af(k, z, 0) and B ± (k } z, 0) are the Fourier coefficients for Blake's no-slip case, allows us to express the Fourier 
components of the velocity field as 



w-^r.A) =w^(r,A)+zw b ^(r,A), 



where 



l-2A-Jwf(r,A) = w a X (r,0), 
l-2A^")w 6 L (r,A) = w£(r,0)- 



(B5) 

(B6a) 
(B6b) 
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Differentiating Eq. (B5) and inverting the Fourier transforms leads to 

(l - 2A w^r, A) = w^r, 0) - 2A w^r, A). 



(B7) 



The solution w (r, 0) is such that V (r) + w (r, 0) is Blake's solution for a no-slip surface, that is, w ± (r, 0) = 
2/i 2 Gf - 2/iGf.f , where Gf represents the source dipole in the z-direction, 



Gf (r) = 



and Gf.f represents a z-dipole of z-Stokcslets [47, 48], 







1 r 



dz \8tt(j, r 3 



d 



G-(r) = — ^ 



G*(r)) . 



Furthermore, if we recognize that 

w 6 ± (r,0) = -2/iGf(r), 
we can integrate each term in Eq. (B7) using Eq. (B6b) and integration by parts to obtain 



w ± (r, A) = - J e- s / 2X [(h + s) Gf - Gf.f ](r + (h + s)e z )d s , 



(B8) 



(B9) 



(BIO) 



(Bll) 



Thus the complete image system for a Stokcslet oriented perpendicular to a partial-slip wall is given by = 
U 1 - + V 1 - + w- 1 , where U 1 - and V 1 - are given by Eqs. (8a) and (11a), and where w J -(r, A) is given by (Bll). Notably, 
this is the same image system as in Blake's no-slip solution; in the partial-slip case, however, the image singularities 
are distributed along a line in the — e z direction, with a magnitude that decays exponentially over 2A. 



APPENDIX C: SOLUTION FOR PARALLEL STOKESLET 

The image system for a Stokcslet oriented parallel to a partial-slip wall is more complicated than for the perpen- 
dicular Stokeslet, but conceptually similar. In this case, Eq. (12) for w'l become w| = and 



wi-X 



dw'i 

dz 



(CI) 



with r\ = x 2 + y 2 + h 2 and a = {1, 2}. In Fourier space, and using Eq. (Dl), these conditions become = and 



wi-X 



~dz~ 



1 ( k 2 {l-hk) + 2k 2 -2*^(1 + MQ \ kh 



Here again we can now solve for wll, using the general formalism given by Eq. (A3), gives 

~k 2 (l-hk)(l + Xk) + 2k 2 (l + 2Xk)] _ kh 



A\ 



2 

~hk~z 

2kxk 2 
hk 3 



(l + Afc)(l + 2Afc) 
(l + Xk)(l + hk) + 2Xk 



(l + Afc)(l + 2Afc) 



-kh 



M = o, 



B 



II 2i&i(hfc-l) 
Wc 2 (l + 2Afc) 



-kh 



These can be written as 



A a {k,Z,\)- 1 + 2Xk + 1 + Xk 



(C2) 

(C3a) 

(C3b) 

(C3c) 
(C3d) 

(C4) 
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and 



where 



q'(k,z,o) 



2kf(hk - 1) 



e- k \ C 2 '(k,z,0) 



2fcifc 2 (/ifc- 1) 

hki ' 



-kh 



Dl(k,z,0) = -^e- kh , Dl(k,z,0) 



4fcifc 2 



-kh 



hk3 ' -av~.~.~> ^3 
As for the perpendicular case, we decompose the velocity field as 

w 11 (r, A) = (r, A) + w|(r, A) + z wjj (r, A) , 



where we now have 



1-2A— )w||(r,A) = wjj(r,0), 



so we get 



(l-2A^)wf(r,A) = wjj(r,0). 
(l - 2A j^) w« (r, A) = w" (r, 0) - A^wjj (r, A) - 2Awf (r, A), 



(C5) 

(C6) 
(C7) 

(C8) 

(C9a) 
(C9b) 
(C9c) 

(CIO) 



where w" (r, 0) is such that V" (r) + w" (r, 0) is Blake's image system for a no-slip surface, that is w" (r, 0) = — 2Gf + 
2/iGf.£ - 2/i 2 Gf . Solving for each of the three terms of Eq. (CIO) using Eqs. (C9b) and (C9c) as well as integration 

by parts allows the flow field to be expressed as u" = U" + V" + w", where U" and V" are given by Eqs. (8b) and 
(lib), and w" is given by 



1 f°° 

w»(r,A) = - J l-G s x +hG s z D x -h 2 G°}(r+(h + s)e z )e-^ds 

f°° r -1 2 

+ 4A / gd (r+(h + s)e z ) e - s / 2A - 1 < 
+ 4A^°°g b (r+(/i + S )e z ) [l - (l + ±) 



ds 



(Cll) 



-s/2A 



where we have defined two new singularities 



gd(r) = gi(r) = ~[w||(r,0)]. 



Their Fourier transforms are given by 



gd(k,z) = 



g&(k,z) = 



,0 e 



k\ k\k2 
87r/i \ k ' k 
(1 - hk) (-k\ -jfeifea 



-kz 



k ' k 



,ik\ ) e 



-kz 



and therefore, by inverse Fourier transforms, we find 



877/x ( Gz r 3 ) 



r~iRD 

-^z;y 



»(r) = (l + ^) [Gf ] = [Gf - fcG« ](r) 



(C12) 

(C13) 
(C14) 

(C15) 
(C16) 
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where Gf.® is the rotlet dipole in the (y; ^-direction, that is, the y-dipole of the z-rotlet Gf , flow field due to a point 
torque and defined as 



C R — -(C s _ c s \ 



and therefore 

r~,RD _ _ (~,SD \ 

z;v 2 v ' xy x >vv'' 
As a summary, the solution for wll in the case of the parallel Stokeslet is given by 

w«(r,A) = l-G s x +hG s z °-h 2 G°}(r+(h + s)e z )e- s / 2X ds 

- 4A / G^(r + (h + s)e z ) 



(C17) 



(C18) 



e s/2X _ 1 



ds 



(C19) 



+ 4A ^ [Gf - hG%]{r + (h + a)e z ) [l - (l + ^ 



['-(" 









ds, 



which is, again, a line integral of fundamental singularities distributed along a line in the — e z direction, with a 
weighted magnitude. Unlike the perpendicular case, however, not all of the magnitudes of the singularities decay 
exponentially away from the image location. 

APPENDIX D: SOME USEFUL FOURIER TRANSFORMS 

Below is a list of two-dimensional Fourier transforms used to derive the image systems in Appendices B and C: 



T 



1 



-kz 



T 

3xy 



- 1^1 P - kz 
~ k 



k 2 kfzk kz 



k 3 



T 
T 
T 
T 



= e 



1_ (:\z 2 
1 



1 



-kz 

i 

- i 

3x 2 



3xz 



ik\ e 



-kz 



= ke 



k\ 



-kz 



-kz 



kik 2 (l + zk) kz 
fc 3 



(Dla) 
(Dlb) 
(Die) 
(Did) 



where we have used the notation r 2 = x 2 + y 2 + z 2 . 



APPENDIX E: FIRST INFLUENCE OF SLIP LENGTH ON PARTICLE DIFFUSIVITIES: 

ALTERNATIVE METHOD 



We show in this section that the results given by Eq. (36) for small slip length can also been obtained by using the 
reciprocal theorem. Let us consider the volume V of fluid above the solid surface S and two steady velocity fields, u 
and u, with corresponding stress tensors, <x and <r, and volume forcing f and f respectively. The reciprocal theorem 
states that 



/ u-ff-e,dS- / u<re z dS= [ fudV- [ fu 
Js Js Jv Jv 



dV. 



(El) 



We take u to be the flow field due to the point force F at (x,y, z) — (0,0, h) near a no-slip surface. Let us also 
consider the flow field, v, due to the same point force near a surface with small slip length, in the sense A h. we 
perform a regular perturbation expansion and write 



v = u + Knu + C(Kn 2 ). 



(E2) 
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The boundary condition for v on the surface is v = Adv/dz, which becomes, at leading order in O (Kn), u = hdu/dz. 
Applying the reciprocal theorem for u and u and using the relation between stress tensor and rate-of-strain tensor 
leads to 



(E3) 



where \i is the shear viscosity of the fluid. In the case of a Stokeslet perpendicular to the surface, F = F z e z , Eq. (E3) 
becomes 



FM Z 



9F? 



2lTflhJ Q (1+u 2 ) 5 



du 



_9Fl_ 



so 



^ = 1-^(1-Kn). 
D Q 8h K ' 

In the case of a Stokeslet parallel to the surface F = F x e x , Eq. (E3) becomes 



F X U X — 



9F? 



A-Kfih J (1 + u 2 ) 5 



du 



9^ 2 
96irnh 



so 



D 



9a 
16ft 



(1 - Kn) ■ 



The results given by Eqs. (E5) and (E7) are the same as those obtained in Eq. (36). 

APPENDIX F: CALCULATION OF COUPLED MOBILITIES 

For convenience, we define below the integrals 

1 f°° 

h = - J l-G s x +hG s z °-h 2 G°}(r+(h + S )e z )e- s / 2X d S , 



I 2 
Is 



4X1 G™(r+(h + s)e z ) 



s/2X _ ^ 



ds, 



= 4A^°°[Gf - hGg ](r + (ft + s)e z ) [l - (l + ±) 



-s/2A 



ds, 



so that the complete image system for a parallel Stokeslet, Eq. (C19), is written wll (r, A) = Ii + I 2 + I3. 



(E4) 



(E5) 



(E6) 



(E7) 



(Fla) 
(Fib) 
(Flc) 



1. Calculation of b X2X1 

Because of the decomposition in Eq. (6), the mobility is given by 

b X2Xl (Kn) = b X2Xl (Kn = 00) + w\ {de x + he z , A). 



(F2) 
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The integrals in Eq. (Fl) and their asymptotic behaviors for small £ are given by 



Ii • e x 



h ■ e x 



I3 • e x 



in^id 
1 

47T/id 

1 



1 r 

Wo 



i/Kn 



C(l + 2u) 



du 



1 



-2 + £ + 5Kn + 4Kn 2 j + O (£ 2 



du 



4:TT(J,d 

1 

iirud 
1 



-^/2Kn + e(Kn+|Kn 2 )+0(^ 



Kn£ 



2 + ^(| — ±u-u 2 ) 



leading to the asymptotic behavior 



[i+£(i + u ) 2 p/ 2 [i + e(i + «) 2 ] 7 / 2 

£ 1/2 Kn - £(3Kn + 4Kn 2 ) + O (£ 2 ) 



du 



wf(de x + he z , A) 



An/id 



-2 + e(§+3Kn+§Kii 9 )+0(0 



and therefore, using Eqs. (42) and (F2), to the mobility 

3£ 



6 X2Xl (Kn) 



(l + 2Kn + Kn 2 ) [1 + 0(0]- 



This asymptotic behavior is valid in the limit d 3> max(/i, A). 



(F3a) 
(F3b) 
(F3c) 
(F3d) 
(F3e) 
(F3f) 

(F4) 

(F5) 



2. Calculation of b v 



In this section, we suppose that the two particles are aligned along the y-axis, at a distance d from each other, 
and apply a force in the x-direction on the first particle. The mobility obtained in this case is the component b X2Xl 
for particles aligned along y, which is equal, by symmetry, to the b y2Vl component for particles aligned along x. 
Consequently, the mobility is given by 



*W ( Kn ) = *W ( Kn = 00 ) + w l ( de y + he z, A). 
The integrals in Eq. (Fl) and their asymptotic behaviors for small £ are now given by 

1 



he x = - 



An^d 
1 



1 r 



[1 + £(!+«) 
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I2 • e x 



I3 • e x 



An/id 
1 

Airfid 
1 

A-K\id 
1 



-1 + d\ + \ku + Kn 2 ) - e (jKn 2 + \n 3 + 9Kn 4 ) O (C< ) 



Kn£ 



Jo 



2-^1 + u) 2 

[l+£(l+ M ) 2 ] 5 /2 



(e~u/Kn_ i y 



du 



( 1/2 Kn - £(2Kn + 3Kn 2 ) + £ 2 ( 2Kn + 9Kn 2 + 21Kn 3 + yKn 4 ) + O (£ 3 



Air/id 
1 

An^id 



Kn£ 



[l + £(l + u) 2 ] 5 / 2 



( 1 -( 1 + i^)^ ,K ") d » 



-£ 1/2 Kn + £ ( ^Kn + 2Kn 2 j - £ 2 f^Kn + 6Kn 2 + yKn 3 + 12Kn 4 ) + O (£ 3 



so that 



w\(de y + he z , A) 



An/id 



+ + e ( §Kn + ^Kn 2 + 3Kn 3 + | K n 4 ) + O (£ 3 ) 
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and therefore, using Eqs. (42) and (F6), the mobility is given by 



b V 2Vi ( Kn ) 



This asymptotic behavior is valid in the limit d » max(/i, A). 



1 + 4Kn + 12Kn + 16Kn J + 8Kn) [1 + O (£)]• 



In this case, we have 



with 



<(r,A) = 



3. Calculation of b X2Z1 (= — b Z2X1 ) 

b X2Zl (Kn) = b X2Zl (Kn = oo) + w^(de x + he z , A), 
?/* [°° 1 + £( U 2 + I U -I) u 



1 

Sir^id 
1 

Sn^id 



[l+£(l+ U ) 2 ] 5 /2 

9, 



/ Kn du 



ill f 

Kn J 

i 1/2 - i V2 ( 3 + ^Kn + 3Kn 2 \ + O (^ 5/2 ) 



and therefore, using Eqs. (42) and (F10), we obtain 

3^3/2 



b X2Zl (Kn) 



lQir^xd 

This asymptotic behavior is valid in the limit d 2 3> max(/i 2 , A 2 ). 



l + 3Kn + 2Kn 2 ) [! + £>(£)]. 



(F9) 



(F10) 

(Fll) 
(F12) 

(F13) 
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